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Abstract 

We study a quantum extension of the spherical p-spin-glass model using 
the imaginary-time replica formalism. We solve the model numerically and 
we discuss two analytical approximation schemes that capture most of the 
features of the solution. The phase diagram and the physical properties of 
the system are determined in two ways: by imposing the usual conditions of 
thermodynamic equilibrium and by using the condition of marginal stability. 
In both cases, the phase diagram consists of two qualitatively different regions. 
If the transition temperature is higher than a critical value T*, quantum ef- 
fects are qualitatively irrelevant and the phase transition is second order, as in 
the classical case. However, when quantum fluctuations depress the transition 
temperature below T*, the transition becomes first order. The susceptibility 
is discontinuous and shows hysteresis across the first order line, a behavior 
reminiscent of that observed in the dipolar Ising spin-glass LiHo x Yi_ x F4 in 
an external transverse magnetic field. We discuss in detail the thermodynam- 
ics and the stationary dynamics of both states. The spectrum of magnetic 
excitations of the equilibrium spin-glass state is gaped, leading to an exponen- 
tially small specific heat at low temperatures. That of the marginally stable 
state is gapless and its specific heat varies linearly with temperature, as gen- 
erally observed in glasses at low temperature. We show that the properties of 
the marginally stable state are closely related to those obtained in studies of 
the real-time dynamics of the system weakly coupled to a quantum thermal 
bath. Finally, we discuss a possible application of our results to the problem 
of polymers in random media. 



1 



1 Introduction 



The description of spin glasses in terms of classical statistical mechanics is generally 
justified since in most cases the transition temperature T g is too high for quantum 
effects to be relevant. In some cases of practical interest, however, quantum fluctu- 
ations, controlled by an external parameter (e.g. magnetic field, doping, pressure), 
may reduce T g down to arbitrarily low values and even suppress the glass transition 
altogether if they are strong enough. Among the experimental systems belonging 
to this class we may cite magnetic systems such as LiHo a; Y 1 _ a .F4[]I[, La 2 _ x Sr :c Cu04 
[0 and UCu5_ x Pd a . H as well as some randomly mixed hydrogen-bonded ferro- 
antiferroelectric crystals The theory of phase transitions in systems such as 
these must necessarily take into account their quantum mechanical nature. 

A question of fundamental interest is whether there exists any qualitative dif- 
ferences between quantum spin-glass systems and their classical counterparts. This 
issue has been extensively investigated experimentally in the case of the com- 
pound LiHo. r Y 1 _ :E F4. This is a site-diluted derivative of the dipolar Ising ferromag- 
net LiHoF 4 . For x < 1, the positional disorder of the magnetic Ho 3+ ions makes 
the long range dipolar couplings random. LiHo x Y!_ x F 4 is thus a spin-glass with 
a freezing temperature T g (x) 0. The application of an external magnetic field H 
transverse to the easy axis allows quantum tunneling through the barrier separating 
the two degenerate ground states of the Ho 3+ ions. Quantum tunneling competes 
against spin freezing and the spin-glass ground state is expected to be destroyed at 
all temperatures if the tunneling frequency (oc H 2 ) is sufficiently high. It was found 
experimentally that LiHo .i67Yo.833F4 is paramagnetic at all temperatures above 
H c (0) « 12 kOe 0. Below this critical field, paramagnetic (pm) and spin-glass 
(sg) phases exist, separated by a line H C (T). Above 25 mK, the phase transition is 
second-order and signaled by a divergence of the non-linear susceptibility ^3- How- 
ever, below 25 mK, the divergence of X3 a t t ne transition becomes a flat maximum 
at the position of which the imaginary part of the low-frequency linear susceptibility 
x"(oj) has a jump. These features strongly suggest that, at low transition temper- 
atures [and in particular at the quantum critical point at (H = H c (0),T = 0)] the 
transverse-field-induced spin-glass transition becomes first order ffl. This conclu- 
sion has recently received further support from the observation of hysteresis in the 
linear susceptibility measured as a function of the transverse field [0. 

From the theoretical point of view, there is no evidence for this type of transition 
in the standard models for quantum spin glasses discussed in the literature. The 
quantum Edwards-Anderson (EA) || and Sherrington-Kirkpatrick (SK) 0, ||, [J 
models are known to undergo second order transitions as also do models of metallic 



spin glasses [10]. Interestingly enough, this type of scenario does occur in less 



standard models with multi-spin interactions [11, O, 13, 14, 15 



In this paper we study in detail the properties of one of these models, the p- 
spin spherical model with random interactions. We solve it using two different 
approaches. The first one consists in imposing the usual conditions of equilibrium 
to find the possible thermodynamic states of the system. In the second one, these are 
determined by imposing a condition of marginal stability. In both cases the problem 
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is solved using a 1-step replica symmetry breaking (rsb) Ansatz shown to be exact 
just as it is in the classical case. We find that, in terms of a quantum parameter T 
to be defined below, there is a line T C (T) in the T — T plane that separates SG and 
PM phases. This line ends at a quantum critical point at T = 0, T = T c (0) above 
which the system is paramagnetic at all temperatures. T thus plays a role similar 
to that of the transverse field in the experiments described above. One of the main 
results of this paper is that a tricritical point (T*, r*) divides the transition line 
in two sections. For T > T*, the SG transition is second order and the behavior 
of the quantum system is in all respects similar to that of the classical one. For 
T < T*, quantum fluctuations drive the transition first order. There is latent heat 
and the magnetic susceptibility is discontinuous and shows hysteresis across the 
first-order line. This is reminiscent of the behavior observed in LiHo a; Y 1 _ 2 ,F4 in a 
transverse magnetic field although, as it will be seen, important differences between 
the predictions of this model and experiment exist. 

It should be noted that other quantum [|TT| , [12j and classical |13|, models with 
similar characteristics were previously discussed in the literature but the connection 
with the experimental results was apparently not realized. 



Several subtle points arise in this and related models [[L2], |13| when one has to 
choose between several possible equilibrium or marginally stable solutions. They 
stem from the existence of multiple stable PM and SG solutions in finite regions of 
phase space. To choose between them one has resort to physical arguments that go 
beyond the usual prescriptions of replica theory. We first show how this problem 
arises within an improved "static approximation" @ that, despite its simplicity, 
contains the essential physics of the model. We then discuss it in the framework of 
the exact numerical solution of the model. 

We investigate the stationary dynamics of the equilibrium and marginally stable 
SG states and show that they are quite different. We find that the magnetic exci- 
tation spectrum of the spin-glass state is gaped, leading to an exponentially small 
specific heat at low temperatures. The spectrum of the marginally stable state is 
gapless and its specific heat varies linearly with temperature. A power-law behav- 



ior of C V (T) is commonly observed in glasses at low-temperatures |16, 17] and is 
explained with models based on a distribution of two-level systems |T^]. In this 
approach there are no such two-level systems but a linear C V (T) stems from the 
fact that the condition of marginal stability selects flat directions in phase space as 
opposed to the equilibrium condition that selects well-defined minima of the free- 
energy. 

We show explicitly that the use of the condition of marginal stability in the 
Matsubara replica approach allows us to obtain, from a purely static calculation, 
partial information about the non- equilibrium real-time dynamics of the same sys- 
tem in contact with an environment, in the long-time, weak-coupling limits (taken 
in this order). We show that the transition line coincides with the dynamic tran- 
sition line determined from the real-time formalism. We find that it also changes 
from second to first order at a tricritical point. This feature was also found in a 
study of the real-time dynamics of the model coupled to a bath where it is signaled 
by a jump in the asymptotic energy density when the system is driven across the 
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transition line [till. We show that m, the break point in Parisi's RSB scheme, co- 



incides with T/T EFF pD |, where T is the temperature of the quantum environment 



(that may be zero) and T EFF is the dynamically generated effective temperature |21 



Finally, we demonstrate that the time-dependent correlation function calculated in 
the replica approach with the condition of marginal stability coincides, in the long 
waiting-time and weak-coupling limits, with the outcome of the dynamic calcula- 
tion for the stationary part of the symmetrized correlation function [p^| . As in the 
classical problem, equilibrium and marginal results can be interpreted in terms of 
the solutions of the TAP equations extended to include quantum fluctuations . 

The paper is organized as follows. In Section [2] we present the model and derive 
its free-energy density with the Matsubara replica formalism. We discuss in detail 
several possible interpretations of the model and its relation to other ones that have 
already been studied. In Section |3] we introduce a refined static approximation which 
enables us to find the correct qualitative behavior of the system in the whole phase 
space. We analyze the consequences of using the static and marginal prescriptions 
for the spin-glass phase within this approximation. In Section |] we present the 
exact numerical solution of the model distinguish again between the equilibrium 
and marginal cases. We also present a low-temperature, low-frequency approximate 
solution, that goes beyond the static approximation, and yields results for the real- 
frequency dependence of equilibrium correlation functions. The numerical results 
are presented in Section [5]. The connection between our results and the real-time 
dynamics of the system is established in Section where we compare the results 
for several quantities obtained in the static and dynamic formalisms. Finally, in 
Section [7] we draw our conclusions and we reinterpret our results for a classical 
polymer in a random media. We also briefly discuss related work in progress. A 
short account of some of our results appeared in Ref . [Ill] . 



2 The model 



We study a quantum extension of the classical spherical p-spin-glass model [£4| in 
which we reinterpret the continuous spins, Sj, as coordinate operators and intro- 
duce canonically conjugate momentum operators, 7Tj. Coordinate and momentum 
operators verify the usual commutation relations 

[si,Sj] = [ir^TTj] = , [^i,Sj] = -iHSij . (2.1) 

The quantum spherical model is then defined by adding to the usual potential energy 
a "kinetic energy" term. The Hamiltonian reads 

7T 2 N 

H[n,s,J) = ^r;+ E Jh-i^-Si,, • (2.2) 
zm h<...<i p 

We denote n 2 = n ■ n, s 2 = s ■ s with 7r = (tti, . . . , tin) and s = (s%, . . . , sn). A 
Lagrange multiplier z enforces the mean spherical constraint 

^M = l (2-3) 

i=l 
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where the angular brackets denote the thermodynamic average. 

The interaction strengths J- n ...i v are taken from a Gaussian distribution with zero 
mean and variance 

{ J h -ip) = 2~/VP r i • ^ 

Hereafter the overline represents the average over disorder. The second term in the 
Hamiltonian is a random Gaussian potential energy with zero average and correla- 
tion 

nm^ = ^y-jf) • (2-5) 

2.1 Replica formalism 

In order to study the static properties of the model we need to compute the disorder 
averaged free-energy density. The replica trick allows us to compute it as 

i i yn _ 1 

Pf = ~ T MZ = -- lim , (2.6) 

where the partition function, Z, is given by 

Z = Tr e~ pH . (2.7) 

In the Matsubara formalism, the disorder averaged replicated partition function can 
be written as a functional integral over replicated periodic functions of imaginary 
time, s (t), with a — 1, ... ,n being a replica index. These functions satisfy s a ((3h) = 
s a (0). In order to decouple the p-interactions in the potential energy we introduce 
in Z n the identity 

1 oc J DQ 5 (NQ ab (r, r') - s a {r) ■ s a (r')) 

DQDX exp — J2 / dr dr'X ab (r, r') (NQ ab (r, r') - s a (r) ■ s b (r')) 
2n ^ Jo Jo 



We shall hereafter use boldface to denote matrices in replica space. The averaged 
replicated partition function can be recast as 



Z n = J DsDXDQ exp (-^eff) , (2. 



with 

% 



1 l ^ fPn , ppn 

~ •S'eff : 



Mr') 



ab 

iN _ rf>* J rf* J , A . . Nf3^ 

+ ^1^J dT J dr Aab{T,T)Q ab (T,T) + —)_^Z a 
ab ® ® o, 

+ ^E/ dr dr'QZiry). (2.9) 
4ft ~T Jo Jo 



5 



We denote with a bullet the usual product: Q*b( T ) — Qabij) • • • Qabij), p times, to 
distinguish it from the operational product. O a b is a short-hand notation for the 
differential operator 

O ab (r - t') = t5 ab 5(r - r') (mJ^ - z)j . (2.10) 

It follows that, at the saddle-point, the expectation value of the order-parameter 
Qab(r, t') is given by 

Qab(r, t') = ^(sa(r)-Sb(r')) , (2.11) 

with Qab( T , T ') periodic in r and t' with period /3K. The mean spherical constraint 
reads Q aa ij, r) = 1 for all r. 

Since we are studying an equilibrium problem, all correlation functions are time- 
translational invariant. They are also symmetric in imaginary time due to the 
time-reversal invariance of the Hamiltonian: 

Qab(r, t') = Q ab (r-r') = Q ab (r'-r) , \ ab (r, r') = X ab (r-r') = X ab (r'-r) . (2.12) 

Using these properties we can simplify the effective action and write it as 

1 1 rPh rf3h 



n 2 1 —r Jo Jo In n 

" " r& n T 2 NB rP n NnB 

Ejf dr\ ab (r)Q ab (r)+ J -^j:l o dr Q*(r) + -fz , (2.13) 



where we have further assumed that z a does not depend on the replica index. In 
the following, we work with the Fourier transforms 



I rpn 

s t (cu k ) = -== / dre WkT s t (r 
\ on Jo 



TMJo 

Si(r) = ^E^ fcr ^(^), (2-14) 



with the Matsubara frequencies given by 



2nk 

w fc = — fc = 0,±l,.... (2.15) 



This implies 

r/3h 



rpn 

Qabi^k) = / dr exp(iu} k r) Q ab (r) , (2.16) 
Jo 

Qab(r) = (phy 1 ^2exp(-iuj k T)Q a b(ujk) ■ (2.17) 
k 
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In terms of the Fourier transformed variables the effective action reads 



k ab 

iN 



% ~ 1 ~ 

-O ab (u k ) + -X ab {u} k ) 
h h 



~ . . Nnp 
s b (uj k ) H —z 



J 2 N{3 



+ 2^ X X A ab (w fc )Q a6 (w fc ) + 

ab ab 



i3h 



E / n dr y—^2exp(-iuj k T)Q ab (uj k ) 



(2.18) 



where we used Q a b{u k ) = Q a b{.~ UJ k)- The functional integration over the functions 
s(uj k ) is quadratic and can be explicitly performed. This amounts to replace the 
quadratic term in the action by 



N 



^Trln[i/3(6(o; A ) + A(a; fc 



(2.19) 



where we took into account a factor (f3h)~ n / 2 that comes from the change in variables 
s a (r) — > s a (uj k ) in the partition function. The trace is to be taken over replica 
indices. 

The effective action is now proportional to N and, in the large iV limit, this allows 
us to evaluate the replicated partition function by the steepest descent method. The 
saddle-point equation with respect to \ a b(u k ) reads 



Q = (O + A) . 



By replacing this value of X ab (uj k ) in ( |2.18j ) one obtains 



(2.20) 



j- 'S'eff 

n 



y E Tr In ((^)- 1 Q) + f E " ¥ E OabMQabM 



ab 



J 2 N(3 



Ah 



ab 



0h 



E/ dr l—J2 ex P(-iu k T)Q ab (uj k )\ + 



Nn(3 



z(2.21) 



hp k 

Finally, the averaged replicated partition function as a function of Q becomes 

= exp (-nNG ) (2.22) 

where 

1 



2Gr 



n 



ETrln((^)- 1 Q) - E ( 1 -^E^Wfe 

k k \ ntl ab 



[UJk 



J 2 f3 rP % ( 1 

E /„ dr \^—J2exp(-iuj k T)Q ab (uj k ) 



2hn 



ab 



(3z 



and the free-energy per spin is 



Pf = \imG . 



(2.23) 



(2.24) 
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The saddle-point equation with respect to the order parameter Q a b{^k) reads 
{Mul + z)5 ab = (CT 1 )^ K) + —£ J rfrexp(^ fc r)Q:r 1 (r), (2.25) 



that transforming back to imaginary time becomes 



- \ - z^j 5 ab 5(r) = Q-J(r) + ^QT\r) ■ (2.26) 

Equation ( |2.26|) together with the spherical constraint 

Qa«(0) = l, -^E ( 5-W = l, (2.27) 

are the equations that characterize the different phases of the model. 

From here on, we shall work with dimensionless quantities. We take J as the 
unit of energy and h/J&s the unit of time. Hence, we redefine the imaginary time 
and Matsubara frequencies as 

r = — => cu k = ^. 2.28 
n J 

The Lagrange multiplier is now given by z = z/J and, consequently, Eq. ( 2.26Q 
becomes 



18P_ 

TP 



+ A Qa b (r) = 6 ab (r) + £ J* dr' V -Q' a r\r - r')Q cb (r') 



(2.29) 



where T = h 2 /(MJ) and we have eliminated all hats in order to simplify the nota- 
tion. 

The matrix elements Q a b{ T ) are the order parameters of the model. The diagonal 
components Qaair) = N^ 1 (s a (r) ■ s a (0)) in Eq. ( |2.29| ) have no classical analog. 
Indeed, in the classical limit < r < %(3 and lim^o Qaa(T~) = Q aa (0) = 1 on 
account of the constraint (|2~lf ). In the quantum mechanical case Qaair) encodes 
information about the equilibrium dynamics of the system through its connection 
to the dynamic local susceptibility x(w), 

. . r°° duo ,,, . exp(— uj\t\) 

Qaa(r) = / -/MHVfll ' 2 ' 30 

J-oo Ti 1 — exp(— puj) 

where x"(a>) = ImxM- It follows that Q aa {i^k) = xML + i e -^ and Q™ = 
C(t)\ t ^_ iT , where C is the real-time autocorrelation function j22|. We shall hereafter 
use the notation Q aa (r) = q^r). 

The off-diagonal elements Q ab are r- independent. This is a general and important 
property of quantum disordered systems in equilibrium proven by Bray and Moore 
in their pioneering work on quantum spin-glasses [|7| . The argument goes as follows. 
The overlap matrix for a ^ b is 



NQa^r.r') = (s a (r) ■ s b (r')) . (2.31) 



8 



Before performing the average over disorder, replicas are decoupled and the thermal 
average factorizes. But the averages (s a (r)) are time-independent for any fixed 
configuration of disorder and so is the average of their product. Consequently, the 
Fourier transform of Qa^bir) = q a ^b is given by 

q a ^b = q a ^bft(3Su k ,o- (2.32) 

We conclude that the RSB solution is confined to the zero mode of the Matsubara 
frequencies and that q a ^b is the quantum analog of the classical order-parameter 
matrix. 



2.2 Relationships to other models 



Besides being a quantum extension of the spherical p spin-glass this model can be in- 
terpreted in other ways. First, if we identify s— (si, s 2 , . . . , sn) as a position vector 
in iV-dimensional space, Hamiltonian (|2.2j ) with the constraint ( [2.3|) describes the 
motion of a quantum particle of mass M constrained to move on the iV-hypersphere 
of radius \/N in the presence of a random potential V(s) with p-dependent correla- 
tions given by Eq. fl2.5| ). 

Second, the well known connection between quantum mechanics in D dimensions 
and statistical mechanics in D + 1 dimensions allows us to interpret this model as 
representing a classical closed polymer in a quenched random medium in infinite 
transverse dimensions [27, [25], [26|. The imaginary time r represents the internal 
coordinate of the polymer of contour length L = (3h. The low temperature limit 
of the quantum model corresponds then to the long length limit of the polymer. 
The N components of the field s denote the N transverse coordinates on the spher- 
ical embedding space. The parameter T is a measure of the linear elasticity of the 
polymer since the first term on the left-hand-side of Eq. Q2.29| ) may be thought 
of as deriving from an elastic energy. Finally, the last interaction term is a con- 
sequence of the average over the quenched disorder. The time-dependent replica 
matrix Q a b{ T ) quantifies the correlation between the position of different monomers, 
r being the distance between them, measured following the polymer direction. The 
replica indices represent different states of the polymer. 

There are also several quantum models related to this one. When p = 2 one 
recovers the model studied in Ref. PSI as well as the disordered quantum-rotor model 
29fl . Quantum Ising spin models with p-spin interactions in a transverse field 

A different quantization 



of Ref. 



have been extensively studied in the literature ITT], O, |30 



rule for spherical spins was proposed in Ref. |3TJ and it was later studied in detail, 
including p-spin interactions, in Ref. fl2|]. The difference in the quantization rule 
results in a different imaginary-time dynamics as the equation of motion contains 
first-order imaginary-time derivatives instead of second-order ones as we have here. 

Notice that if we set p — 4, Eq. ( |2.29| ) becomes similar to Eqs. (4)- (6) in Ref. [0 
for the replicated boson Green's function G a b{r) in the SU(A/") Heisenberg spin-glass 
model. Indeed, the equations look identical if we formally take 5* = — 1, redefine 
J and replace the term linear in the Matsubara frequencies by a quadratic one, 
i.e., if we use first-order instead of second-order imaginary-time derivatives. Note, 
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however, that the symmetry properties of the quantities G (in Ref. |32[]) and Q 
(here) are quite different. 

Chandra et al. showed the equivalence between a model for a classical array of 
Josephson junctions in a magnetic field and the p-spin-glass model with p = 4 [38]. 
Using a diagrammatic expansion technique similar to that used to derive classical 
or quantum TAP equations, Kagan et al. |37j showed that the equations governing 
the equilibrium behavior of the quantized model are similar, though not identical, 
to the Matsubara equations we analyze here. 



3 The static approximation 



Before presenting the full numerical solution of the model it is useful to study it 
using a very simple approximation known as the static approximation. The latter, 
first introduced as a variational Ansatz in Ref. ||, consists in proposing qdij) = q d , 
independent of r, and determining by minimization of the free-energy. Despite its 
simplicity, detailed comparison with exact numerical results || [12| has shown that 
the static approximation captures much of the physics of quantum spin mean-field 
models. In our case, however, the naive static approximation applied at the level 
of the saddle point equations is oversimplified since the spherical constraint ( |2.27| ) 
immediately sets q^ = 1 thus loosing all variational freedom. We then proceed 
differently and replace qdij) by qd and \dij) by at the level of the effective action 
(|2.9|) and treat both quantities as variational parameters. In this way, we do not 
impose the spherical constraint strictly (since qd can be different from 1) but in an 
averaged manner. 

The effective action in the static approximation becomes 



^ ^EFF 

h 



N 



TrlnQ 



Nn 



£1* 







UJt 



+ z 



ab 



Nn 



(3.1) 



By using the identity |33 



£1* 

k 




2 In 



2 sinh 



we get the following expression for the free-energy density 



Pf 



— \n((3z) + In 
2 



2 sinh 



— lim — 

n->0 fl 



--(l + /3z(l-g d )) . 



The saddle-point equation with respect to z reads 



(3.2) 



iTrlnQ + ^g: 

ab 



(3.3) 



1 + -J- - -\ - coth ( -VTz 
f3z 2V z \2 



(3.4) 
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and replaces the spherical constraint. As expected, it is independent of the non- 
diagonal part of the Q matrix. At high temperature (/3 — > 0) one recovers the 
classical result, q d = 1. In the opposite limit, /3 — > oo, 

(3.5) 

The second term on the right hand side of the equation corresponds to the zero-point 
reduction of the order parameter characteristic of quantum magnets. 
The saddle-point equation with respect to the order-parameter Q ab is 

= \q-j + y^r 1 • (3.6) 

3.1 The paramagnetic solution 

At high temperatures or if quantum fluctuations are strong we expect a PM phase. 
Within the static approximation, the PM solution is associated to a diagonal replica 
matrix Q: 

Qab(r) = q d 5 ab . (3.7) 
The saddle-point equation ( |3.6| ) is now very simple 

1 I PP P-I (O fi\ 

» = Wd + T ti • (3.8) 

This equation and Eq. (|3.4j) form a set of two coupled equations for q d and z that 
can be solved numerically for all f3 and V. 

Some asymptotic behaviors can be obtained analytically. First, in the classical 
limit, T fixed and r — > 0, Eqs. ( |3.6| ) and ( |3.8| ) imply q d = 1 and (3z = 1 +p/3 2 /2, 
in agreement with the results in Ref. [|4jj. Second, when T ^> 1 but f3Y <C 1, 
quantum fluctuations are still irrelevant and we get (3z = 1 and q d = 1. Third, in 
the limit r — > oo and T finite, q d tends to zero as q d ~ 4/(/3r) and z ~ T/4. The 
zero-temperature limit is more subtle and we discuss it below. 

The properties of the general solution of Eqs. ( |3.6| ) and ( p.8[) are best understood 
graphically. In Figs. [l|(a) and (b) we display g(q d ) = q d — rhs(Eq. ( p.4|) ) vs. g^, 
obtained replacing 2 by its value given by Eq. ( |3.8| ), for several values of T and T. 
The points where g(q d ) crosses the horizontal axis are the solutions for q d . 

The summary of our results for p > 3 is the following (the case p = 2 is different 
and we discuss it in Section |3~3"1) : 

• For fixed f3 < /3 P (« 6 for p = 3) and all values of T, the equations admit only 
one solution with < q d < 1. The function q d (T) decreases monotonically 
with q d (0) = 1 and q d — > when r — > 00. This is represented in Fig. |l| (a), 
where we show g(q d ) for /3 = 1 < @ P and T = 0.1, 10, 100. The analysis of the 
free-energy density shows that this solution is a local minimum of f(qd)- 
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_Q g ill I ■ I i I i I i I i 

'0.0 0.2 0.4 0.6 0.8 1.0 

Figure 1: (a): The function g(qd) that determines the PM solutions for 
P = 1 < fa « 6 and T = 0.1,10,100, from bottom to top. For all T we 
find only one solution that moves towards smaller values of qj when V increases, 
(b): The same for (3 = 10 > f3 P ~ 6 and several choices of T. For T = 6 > r cl (top 
curve) there is only one solution with a small value of qa, for T c2 < T = 4 < T cl 
(middle curve) there are three solutions, while for T = 2 < T c2 (lower curve) there 
is again only one solution but with a large value of qj. 
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• For (3 > /3 P , there are several PM solutions and phase transitions can occur 
between them upon varying T. For T above a critical value r cl , there exists 
a solution PMi with a small value of qd = q^ (the upper curve in Fig. |l](b) ) 
that we call a quantum paramagnet. This solution can be continued below r cl 
but two new ones appear: PM 2 , that we call a classical paramagnet, with large 
Qd — Qd i an d PM 3 with an intermediate value of qd- The three solutions coexist 
in the interval r c2 < T < F c i- At r c2 , PMi and PM3 merge and disappear. The 
three coexisting solutions are shown in the intermediate curve in Fig. |l](b). 
Below r c2 , only PM 2 exists (the lower curve in Fig. 0(b)). An analysis of the 
stability of these solutions shows that, where they exist, PMx and PM 2 are local 
minima of the free-energy whereas PM 3 is a local maximum and can therefore 
be immediately discarded. In the limit (3 — > 00 the two stable solutions can be 
simply computed and one obtains q^ ~ 1 — yT 1 / (2p)(3~ 1 / 2 — > 1 with z ~ p/3/2 
and q^ ~ 4/(T/3) — > with z ~ T/4. Both solutions exist for all values of T 
when T — > 0, indicating that r c i — > 00 and r c2 — > as T — > 0. The dotted 
lines in Fig. |5] delimit the region where there is coexistence of paramagnetic 
solutions in the case p = 3. 

• The free-energies of the two stable PM solutions cross at r c0 , a coupling inter- 
mediate between T cl and T c2 . For T > T c0 , f(qf) < /(q^) but the inequality 
is reversed for T < T c0 . This would suggest that, for all < T < T P , there is 
a first-order transition between the two stable PM states with PM 2 favored for 
T < r c o. However, we will see below that this transition is only possible for 
T* < T < T P where T* is slightly below T P for small to moderate p. 



3.2 The SG solution 

3.2.1 The replica symmetric solution 

At low temperatures one expects the system to exhibit a non trivial SG solution 
characterized by a replica matrix with non-vanishing off-diagonal elements. The 
simplest possible Ansatz is the replica-symmetric (rs) one, 

Qab = (qd - <?ea)<U + <?ea- (3.9) 



From the general expression ( |3.3|) for the free-energy one then derives three 
equations for z, qd and q EA . The equation for z is independent of the Ansatz and 
hence coincides with Eq. Q3.4j) . The other two equations read 

1 q d - 2q BA (3p ! 

+ ^r<fd , I 3 - 10 ) 



(3 (q d - g EA ) 2 2 
1 , p(3 



2 



after we eliminated the paramagnetic solution q EA = 0. 

These equations have non-trivial solutions in a finite region of phase space satis- 
fying < q EA < qd- The phenomenon of multiplicity of solutions at low temperatures 
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found in the PM case appears also here. For temperatures that are not too low, (e.g. 
for P — 10 for p — 3), these equations have two solutions, one stable and one unsta- 
ble within the RS Ansatz. At very low temperatures, two other solutions appear, one 
stable and one unstable within RS. However, a careful analysis of stability shows 
that, as in the classical case the RS solutions are always unstable with respect to 



replica symmetry breaking. To show this, it is sufficient to compute the replicon At, 
the transverse eigenvalue of the Hessian matrix . We will derive a general expression 
for At in Section |4.4.1|. Here we simply quote the result in the static approximation 



which is 



At = - - ^P(P - 1K~ A 2 . (3.12) 



It follows from this equation and Eq. Q3.1ip that cannot be positive for any of 



the RS solutions if p > 3. (See Section |3]3| for a discussion of the case p = 2.) 

3.2.2 The 1-step replica symmetry breaking solution 

Using Parisi's RSB scheme, we search now for a 1-step (rsb) solutions of the form 

Qab = (qd - <?e A )<U + (<? EA - goKft + go , (3.13) 

where the matrix e is defined as 

J 1 if a and b are in a diagonal block , , 

6ab ~ { otherwise [ } 

We show in the Appendix that this Ansatz actually gives the exact solution to the full 
problem. In the absence of an external magnetic field go — and the extremization 
equations with respect to qd, g E A arid m read 

q d + (m - 2)q EA p(3 2 p_j 

P z = 7 r-i yu \ + ^~ q d ' 3 - 15 

= t TT^-Wf V-^fc" 2 , (3-16) 

(qd - 5ea(1 - m))(qd - g E Aj 2 

n Pi/' qd-q^A \ . Pf3 2 „ 

= —7 t: rr + — In r + — C , (3.17) 

m(q d - g BA (l - m)) \g d - g EA (l - m) J 2 



where we excluded the two solutions m = 1 or g EA = from Eq. ( |3.16| ) . The first and 



second equations above have new m-dependent factors. These equations reduce to 
the RS equations when m = 0. The third equation is new and is used to determine 
the break point m. 

It is useful to define the parameters 

y = — , x p = -^-, (3.18) 
qd i-y 
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Subtracting Eq. ( |3.17| ) from Eq. (|3.16| ) and writing the result in terms of x p one 
obtains 




, . =0. (3.19) 

1 + x p p(l + Xp) 

This is a master equation for x p that only depends on the parameter p. Once x p is 
known from this equation, y follows from y = x p /(m + x p ). For p = 3, a case that 
we will consider in detail 

x 3 = 1.81696 (3.20) 

Equation ( p,16|) yields a simple relation between q EA and x p . Using then q EA = 
Xp/ (m + x p )qd one has 

p^)\ v _ _i_ (321) 

P(m(3) 2 xl->(m + x p y 



For fixed values of m and /3, one obtains q EA , q^ and z from Eqs. fl3.21|) , (|3.22j ) 



and ( |3.15| ), respectively. Solving numerically Eq. ( |3.4| ) one gets the corresponding 



value of T. 

To study the transitions between the SG and PM phases described above we 
compare their free energies. Three cases must be distinguished 

• (3 < /3 P (but greater than the inverse of the classical critical temperature). 
This case is represented in Fig. ^] for p = 3 and j3 = 4. The SG and PM free 
energies increase monotonically with T. At T = T c , f SG = f PM and a phase 
transition occurs. 

As in the classical case (where the temperature is the control parameter), 
m = 1 at the transition. The spin-glass solution continues to exist beyond the 
transition point until T = T MAX > T c where it disappears. However, between 
T c and r MAX , m > 1 . Values of m greater than 1 do not represent physically 
allowed states and should not be considered. This follows from the fact that 
the susceptibility x — ft [Qd ~ (1 — m )lEA] can be shown to be strictly smaller 
than f3qdi which implies < m < 1. Below the critical point, / SG > / PM 
meaning that the SG solution maximizes the free-energy. 

The PM solution exists below the transition all the way down to T = 0. This 
solution is locally stable for all values of T and could be interpreted as a 
metastable state below T c . We believe, however, that this solution should 
be discarded on physical grounds because it is continuosuly connected to a 
T = state whose thermodynamic properties at low temperatures are un- 
physical: the internal energy and the susceptibility of this state diverge as 
T — > whereas these quantities can be shown to be finite in the ground state 
of Hamiltonian ( |2.2| ). Therefore, we shall consider that the continuation of 
the PM solution below T c is a spurious solution and that the only physically 
allowed phase below r c in this temperature range is the SG. 
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The Edwards- Anderson order parameter g EA is discontinuous at the transition. 
The transition is nevertheless of second order in the thermodynamic sense. 
This follows from the fact that m = 1 at r c and so the effective number of 
degrees of freedom involved is (1 — m)q EA — > at r c . Therefore, there is no 
latent heat and the linear susceptibility is continuous. However, in contrast 
with the SK-model, the transition is not associated with a divergence of the 
SG susceptibility. The situation is the same as in the replica theory of the 
classical model [ fMf . 

ft* < ft < ftp where T P — T* <^ T P except for very large values of p. The case 
p = 10 for ft = 6.5 is represented in Fig. |3] where we show the free-energies 
of the two PM and the SG solutions as functions of T. Decreasing T from 
r > 1 we encounter two phase transitions. The first one occurs when the free 
energies of the two PM phases cross which signals a first order phase transition 
between the quantum and classical paramagnets. The glass transition occurs 
at a lower coupling, when the free energies of the PM2 and SG states cross. 
This transition is similar to that found for T > T P and it is characterized by a 
discontinuity of the order parameter and the value m = 1 at the transition. As 
in the previous case, / SG > f PM2 and the SG solution maximizes the free-energy. 

ft > ft*. This case is represented in Fig. |] for p = 3 and ft = 10. The free- 
energies of the states PMi and PM 2 cross as before. However, / SG does not cross 
f PM2 but it crosses / PMl instead. Therefore, a transition to the low-temperature 
ordered state can occur from PM 1; not from PM 2 . If we were to interpret the 
crossings of the free energies in the standard way, we would conclude that below 
the transition to PM 2 the system stays in that state down to T = 0. However, 
PM 2 , which is the continuation to low temperatures of the PM state discussed 
above for T < r c and T > T P , exhibits the same unphysical properties. We 
thus conclude that, for ft > ft*, the PMi to PM 2 transition is spurious and that 
the system stays in the PMi state until the latter disappears in favor of the SG 
state. It should be noted that a similar situation was encountered in a study of 
the Ising SG model in a transverse field at low temperatures [31|. In contrast 
to the situation for ft < ft*, f SG < f PMl and the SG solution minimizes the 
free-energy Such an inversion of the order of the free-energies was also found 
in a study of an anisotropic classical p-spin model JJS). It can be thought as a 
result of the competition between the conflicting requirement of maximization 
and minimization with respect to q EA and q^, respectively. 

We find that the break point m considered as a function of V has two branches. 
Physical values of m must be chosen from the branch that can be continued to 
the classical limit, T = 0. For ft > ft*, the maximum physical value of the break 
point m MAX < 1. The SG solution ceases to exist at the corresponding value of 
r, r MAX (this will become more transparent in the zero temperature limit that 
we discuss next). An important consequence is that now m is discontinuous at 
the transition implying that the latter becomes first order. Since the physical 
SG and PM solutions extend beyond the point where their free-energies cross, 
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Figure 2: The SG and PM free-energies (solid and dashed lines, respectively) in the 
static approximation as a function of T for (5 — 4 and p = 3. 

there is a region of phase coexistence and hysteresis is expected (see below). 

The phase diagram resulting from this analysis is represented for p = 3 in Fig. [5] 
(thin lines). The PMx to PM 2 transition is invisible on the scale of the plot for this 
value of p. At low temperatures, the first order line exhibits reentrant behavior, as 



is the case for the Ising p-spin-glass model in a transverse field ||12|| . This feature 
is an artifact of the static approximation that disappears in the exact treatment of 
the problem, as we shall see in Section f|. The second and first order critical lines 
meet at the tricritical point (T*,T*). It is interesting to notice that the 180° rule 



36 that imposes that the angle formed by the intersection of two critical lines at a 



tricritical point should be < 180° is satisfied in the static approximation. 

3.2.3 The static approximation in the zero temperature limit 

The behavior of the system can be completely elucidated at zero temperature. In 
this limit, the requirement that q EA is finite implies 

(3m < +oo , (3.23) 

meaning that m — > 0. This implies that Eq. ( |3.22j ) can be simplified to 



~ ^ (3.24) 



and coth(/9\/r2;/2) in Eq. ( |3.4| ) can be replaced by one. The Lagrange multiplier in 
Eq. (|3.15| ) becomes 

z ~ X \ - (3.25) 

(1 + x p )q d f3m 
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0.5 




Figure 3: The SG free-energy (solid line) and the free-energies of PMi and PM 2 (long- 
dashed and dashed lines, respectively) in the static approximation as a function of 
r for (3 = 6.5 and p = 10. In this case / SG intersects f PU2 . 




2 4 6 

r 

Figure 4: The SG free-energy (solid line) and the free-energies of PMi and PM2 (long- 
dashed and dashed lines, respectively) in the static approximation as a function of 
r for P = 10 and p = 3. In this case / SG intersects / PMl . 
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Figure 5: Static (thin lines) and dynamic (thick lines) phase diagrams of the p-spin 
model for p = 3 in the static approximation. Solid and dashed lines represent second 
and first order transitions, respectively. The dotted line delimits the region where 
two stable paramagnetic solutions coexist. 



and we obtain for T as a function of qd, 



8pX% (p_ 2 )/2 /, n2 

—Qd 



\ l + x p 



(3.26) 



The right hand side (rhs) of this equation is a bell-shaped curve, that vanishes at 
qd = and qa — 1. It reaches its maximum at g^ IAX = {p — 2)/(p + 2), independently 
of x p . The corresponding value of Y is 



ipXp 



p-2 
\J 1 + x p Vp + 2 / 



(p-2)/2 



J> + 2 , 



(3.27) 



The physically meaningful solution for T < must be searched on the branch 

corresponding to q d > q^ iAx . This is the branch that gives the correct classical limit, 
qd — 1 for T = T = and has the expected property that q EA is a decreasing 
function of T. This solution can be continued until T = T^™™. However, the 
transition occurs at a lower value r c , the point at which the free-energies of the SG 
and PMi states cross. These are respectively given by 



fsc — 



1 



2 pm 



\n(l + x p ) + 



Tz z , 
2" + 2^ 



1) 



P(3m p 



and 



PMl 



^ (3.28) 



(3.29) 



At r c , q EA is finite and m — 0. The transition is thus of first order. 
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3.2.4 The condition of marginal stability in the static approximation 

The condition of marginal stability leads to a different equation for the break point 
m. In this prescription one does not require m to be an extremum of the free- 
energy; instead, one requires that the replicon eigenvalue vanishes. In the static 
approximation, the expression derived in Section [4.4| yields: 



1 (3 



2 



^(p-l)C- 2 = 0. (3.30) 

The other extremal conditions remain unchanged. Combining Eqs. ( [3.16|) and (|3.30|) , 
and using the definitions ( |3.18| ) we arrive at 

1 - V 

m=(p-2) x p = p-2. (3.31) 

y 

For reasons that will become apparent in Section ||], we define the dynamic tran- 
sition line as the boundary of the region in the T — T plane, with T < T*, where 
the marginally stable SG exists. The line r^(T) consists of two sections. The first 
one, that starts at the classical critical point, (T<j, 0), is determined by the condition 
m = 1. The dynamic transition in this case is of the same nature as in the classical 
case. The condition m = 1 can be fulfilled until reaching a tricritical point (T^, T^). 
At lower temperatures the marginal solution disappears before reaching the value 
m — 1. 

The dynamic phase diagram for p = 3 is shown in Fig.|5| (thick lines). As in the 
equilibrium case, m is discontinuous across the dashed line. Yd lies always above T c . 
According to the interpretation of this line given in Section |6|, this means that the 
equilibrium state can never be reached dynamically starting from an initial state in 
the PM phase. The two lines come very close to each other for T ~ T*. Within the 
accuracy of our calculations we cannot assert whether they precisely touch at T*, an 
intriguing possibility. For T < T*, m varies continuously along T^(T) and vanishes 
at the quantum critical point. 

It is interesting to note that the 180° rule does not hold in this out of 
equilibrium situation. Indeed, in Fig |5] one sees that the first order and second order 
line join continuously at the tricritical point. 

3.2.5 The marginal stable solution in the zero temperature limit 

Again in this case it is simple to obtain an analytic equation for ^(r) at zero 
temperature. Inserting x p = p — 2 in Eq. ( p.26|) leads to 



r = ^8p^-2) !(1 _ ?(i)2?r)/2 (332) 

The maximum of this curve is also located at g^ AX = (p— 2) / (p+2) and the physically 
meaningful solution is on the branch to the right of the maximum. This equation 
has a solution up to 
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At zero temperature, the 1-step RSB solution with the marginal condition disappears 



at T = rj^ G where m is zero. 

Note that if we compare the maximum value of V for which there is a 1-step RSB 
solution with the equilibrium and marginal criteria we have 

pSTATIC v. pMARG jn 

1 MAX ^ 1 MAX • {O.O^i) 

Hence the marginal solutions disappear at lower values of T than the equilibrium 
solutions. This result is also obtained with the quantum TAP approach |23) when 
one follows metastable states as functions of their energy density. This behavior is 
similar to that of the classical problem with varying temperature: the equilibrium 



solution disappears at higher temperatures than the marginal solutions This 
does not mean that the static transition occurs at a higher value of T at fixed 
temperature. In fact, the static transition is determined by the identity of the free- 
energies / SG = / PMl that at zero temperature yields r™ TIC ~ 3.61 and it is smaller 
than the dynamical critical value T° YN = r^| G — 3.97. 

The inequality ( |3.34| ) is valid at finite temperatures and beyond the static ap- 
proximation, too. This can be checked via the full numerical solution of the problem 
or with the quantum tap approach used in Ref. 

3.3 The simpler p = 2 case 

The model with p = 2 has been studied by several authors [2~S, [SPj and already at 
the classical limit it is known to show important differences with respect to the case 
p > 3. The transition is always second order and it is accompanied by a divergence of 
the SG susceptibility. There is no multiplicity of fully stable paramagnetic solutions 
and the spin-glass phase is RS. This statement follows from the fact that, for p = 2, 
Eq. ( p,19| ) only admits x 2 = as a solution. Then, Eq. ( |3.18| ) implies that, for 
q EA 7^ 0, m = 0. Notice that, if p = 2, the replicon eigenvalue (|3.12| ) is identically 
zero. The RS solution is thus marginally stable, as in the classical case 



4 The exact solution 

The static approximation neglects the imaginary-time dependence of the diagonal 
elements of the order-parameter matrix. In this Section we discuss the exact equa- 
tions for the paramagnetic and spin-glass phases. We study the properties of the 
SG phase using both equilibrium and marginality conditions. In the latter case, we 
derive exact equations for the Edwards- Anderson parameter and the break point m 
that we shall later compare to the results of the real-time calculation (Section |6]). 
We present a low-temperature and low-frequency approximation that allows us to 
show that the marginally stable state is the only gapless solution. It also allows us 
to compute the dependence of gd(0) on T, the T-dependence of the specific heat, 
and the low-temperature properties of — ► 0). 
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4.1 Equations for the exact paramagnetic solution 

In the paramagnetic phase the order parameter matrix is diagonal, 

Qab{u k ) = q d (uj k )5 ab . (4.1) 
Replacing in Eq. (|2.25|) we obtain the following equation for qd(u)k) 



f + *= + , (4-2) 

where we have defined the self-energy 

S(wjfe) = | / dr exp(iu k r) gJpV) . (4.3) 
2 Jo 

The Lagrange multiplier z is determined by the spherical constraint 

gd(o) = i£&M = i. (4.4) 

^ k 

Thus, in the general case, z is determined by a non-trivial implicit equation in 
contrast with its analog in the static approximation, Eq. ( |3.8| ). It follows from 
Eq. ( |4.2| ) that the Fourier components of ?d(r) with k ^ are strongly suppressed 
in the limit r — > 0. The normalization condition then implies that limr^o q~d{^k) = 
f35u kt o which is the classical result. 

4.2 Equations for the exact spin-glass solution 

The RS Ansatz in Fourier space is given by 

Qabi^k) = - <?ea) S ab + q EA . (4.5) 

The inverse matrix, Q -1 , can be easily evaluated to yield (in the limit n — > 0) 

(q- 1 ) ( Uk ) = ^-2 Sab ^— 2- (4-6) 

v /ab 9dK)-&A (q d (uJk) - q EA ) 



Inserting the above expressions in Eq. (|2.25| ) and using the identity (|2.32|) , we get 
(for a ^ b and u k = 0) 



— 2 + = • (4.7) 

(&(0) - /%a) 2 2 yEA 1 J 

Note that in the classical limit q~d(0) = (3, and one recovers the classical equation for 
q EA in the RS case ||24j| . As in the static approximation, the replica symmetric Ansatz 



is unstable towards replica symmetry breaking. The replicon eigenvalue associated 
with Eq. O is 

J? (3 2 

(q d (0) - f3q E 



At = t=t J Q _ V2 ~ ^P(P ~ 1)C" 2 • (4-8) 



iEA ; 
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This equation and Eq. ( |4.7| ) yield A T = p{2 — p)q^ 2 ^ /2. This eigenvalue is negative 
for all p > 3. This excludes a RS stable SG state and the symmetry between the 
replicas must be broken. 



Inspired by the classical case [^4| , and the results from the static approximation, 
we use a 1-step rsb Ansatz that we prove to be exact in the Appendix. In Fourier 
space this can be written in the form 



Qabi^k) = (Qdi^k) ~ <?ea) $ab + (<?EA - <?o) tab + % 



(4.9) 



where the matrix e has been defined in Eq. ( |3.14| ). In the absence of an external 
field, the saddle point equation for q yields q = 0, as in the static approximation 

. The inverse matrix, Q _1 (u;fc) can then be written as 



and the classical case 



A(to k )6 ab + B(u k )e 



lib 



where (in the limit n — > 0) 



A(u k ) 



(4.10) 



(4.11) 



and 



B(u k ) 



QdM - qL( m - 1) + qd{^k)q^{m - 2) 

The saddle point equations for q EA and qaiyJk) can be obtained by inserting Eqs 
and ( |4.10|) in Eq. Q2.25| ) at a ^ b, to k = and a = b, respectively. They read 



(4.12) 



Pp-2 



and 



LOT 



+ Z 



«2(0) - P 2 qL(m - 1) + Pq EA q d m™ - 2) 2^ EA 

EA (m — 2)S Wkt Q 
qj(td k ) - P 2 ql A (m - l)6 Ukfi + (3q EA q d {u k )(m - 2)5^ 







(4.13) 



+ t(u k ), (4.14) 



with S(c<jfe) defined in Eq. (|4.3| ). Equations Qj.l3Q and (|4.14|) must be supplemented 
by an equation for the break point parameter, m. As within the static approxima- 
tion, the latter may be obtained using two different prescriptions that we discuss 
next. 



4.3 The equilibrium spin-glass solution 

In equilibrium, the free-energy per spin, /, must be stationary with respect to 
variations in m. The only m-dependent terms in / are 



ili m I('-lTrln(r 1 Q(0))- 



2 



n{m - l)ql A 



(4.15) 
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which can be easily evaluated with the result 



1 

2P 



In 



<7d(0) - g BA (l ~ m) \ m-l ^ 



Qd(0) ~ Qe 







m 



^d(O) - <?ea(1 - m) 



~ -£(™-l)q : 



p 

EA ' 



The extremization equation df /dm = is 

_1 PQeA J_ ^ / gd(0) - /?g E A 

m q d (0) - /3g BA (l - m) m? \ q d (0) - f3q EA (l - m) 



(4.16) 



Combining Eqs. ( |4.13| ) and (|4.17|) and defining 

PQea 



y 



X Ti 



my 



iyC = 0. (4.17) 



(4.18) 



we obtain the same expression for x p that we derived in the static approximation, 
Eq. ( |3.19| ). Using Eq. ( |4.13| ) and the definitions ( |4.18| ), we have 



p(/3m) 2 q 



2 aP 

EA 



2j 

1 ~\~ OCf. 



(4.19) 



and 



1+Xr 



(4.20) 



The first of these equations is identical to Eq. ( |3.21 



It is convenient to separate q d and the self-energy into constant and r-dependent 
parts, 



<M r ) = Qea + 9reg(t) , 
S(r) = ^^ 1 + E REG (r) 



(4.21) 
(4.22) 



Substituting in Eq. (|4.14 ) and using Eqs. ( 4.19|) and (|4.20 ), the terms that are 
proportional to 8^,0 cancel and we obtain an equation for the regular part of q d (u k ), 



l 



<?REG {Uk) 



+ S REG (^fc) — E REG (0) 



where 



V n p-1 1 + X p 



(4.23) 



(4.24) 



and 



S R eg(^) - Shbo(0) = ? / (cos(w fc r) - 1) (gg V) - g^ 1 ) . 

2 jo 



(4.25) 
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4.4 The marginally stable spin-glass solution 



In this Section we compute the replicon eigenvalue for the exact problem and we 
derive the consequences of using the condition of marginal stability that corresponds 
to setting it to zero. 

4.4.1 The replicon eigenvalue 

To derive the replicon eigenvalue, we first compute the second order variation of Go 
with respect to X ab {oj k ) = Q a b(u k )//3 which gives the Gaussian X fluctuations: 

25 2 G = ^TTrfx-^MX^ 



P-2 



1 ab uj' k w'l V ^k 

x exp [-it (u' k + u'i)] 8X ab {u' k )5X ab {ul) . (4.26) 



For the 1-step RSB Ansatz the eigenvalues and eigenvectors of this quadratic form 
may be obtained by solving the equation 



A 2 (uj k )5X ab (uj k ) + A(oo k )B(uj k ) [(SX(w k ) ■ e) ab + (e • 5X(u k ) 



ab 



+B 2 (u k ) (e • 5X(u k ) ■ e) ab - P ^ P 1 \ v- 2 e ab 5X ab (u k ) = (3" 2 A T (u k )5X ab (u k ) (4.27) 



where a ^ b and A(uj k ) and B(uj k ) are defined by Eqs. ( |4.11| ) and (|4.12|) . Since we 



are looking for the replicon (transverse) eigenvalue we impose the constraints 

(e-6X{u k )) gb = 0, (4.28) 
[l-e ab )5X ab (uj k ) = 0, (4.29) 



and set u k = to obtain 



At = /3 2 A 2 (0)-^p( P -1)^- 2 . (4.30) 



The condition of marginal stability corresponds to = and replaces ( [4.17 ) for 
the determination of the break point. 

4.4.2 The block-size m 

Setting A T = and using Eq. ( |4.13| ) we obtain 

m = (p-2) = (p-2) — . (4.31) 

pqea y 



It follows from this expression and Eq. Q4.2C ) that x p = p — 2. This is the result 



that we obtained in the static approximation. In the classical limit, m = (p 
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2)(1 — q EA )/q EA , a known result []24| . Moreover, using again Eq. Q4.13Q , the above 
expression can be written as 



(3 m =(p-2) ] J- ( — T)q -^, (4.32) 



which is identical to Eq. ( |3.21| ) obtained in the static approximation. 



Equations ( |4.23|) - (|4.25|) , derived for the equilibrium state, also hold for the 
marginal SG provided we make the substitution x p — > p — 2. 

4.4.3 The Edwards- Anderson order parameter 



If we rewrite the denominator in the first term of Eq. ( |4.13j ) and replace m by its 



expression in Eq. (|4.32| ) we obtain a new equation for q EA that reads 

P -K^ (UO) ~ /W = ^^C- 2 XL G (0) • (4.33) 



Note the absence of explicit dependence upon T both here and in Eq. ( 4.32 ). The 
quantum parameter only enters implicitly through (jd(0) and q EA . 

An identical equation for the Edwards-Anderson parameter can also be derived 
with the quantum TAP approach |23[ when the threshold value for the free-energy 



density is chosen. In the classical limit, qd(0) = (3 and this equation becomes the 
equation for q EA of the classical TAP states of the threshold level K3|. 



This equation may be studied in the usual way, i.e., by looking at the shape 
of the function that appears on its right-hand-side. This vanishes at q EA = and 
q EA = q~d(0)/P and has a maximum at (3q EA = ^^^(0). The maximum is at m = 2 
just as in the classical problem. Two solutions for /3g EA , one on the right-branch, the 
other one on the left-branch. It is simple to see that the correct physical behavior of 
a decreasing Edwards- Anderson parameter with increasing temperature is achieved 
by the solution on the right-branch. 

Note, however, that other equations where q~d(fy appears need also be satisfied 
in the quantum problem. These depend explicitly on the quantum parameter T. 
Depending on T and T the marginally stable SG may disappear before reaching the 
m = 2 value. 



4.4.4 The low-frequency limit 

The existence of a vanishing transverse eigenvalue implies that the spectrum of 
magnetic excitations of the marginally stable SG state is gapless. Therefore, at 
T = 0, g RE G( r ) is expected to decay asymptotically as a power- law: 

?reg(t) ~ 7—r- |t| > OO . (4.34) 
\ T \ 

Equivalently, the imaginary part of the susceptibility x"{ u ) vanishes as a/ 0-1 ) as 
uj — » 0. In the zero temperature limit the Matsubara frequencies are continuous and 



26 



Eq. ( |4.25| ) becomes 



Sreg(^) - S REG (0) 



p 



dr (coscjt — V 



(p - l)q p E fq REG (r) + 



(4.35) 



The terms represented by the dots contain higher powers of Q'reg (^") ■ I n the us —>■ 
the integral is dominated by the long r limit, and we may replace Qreg (^~) by its 
asymptotic form ( f4.34p . Eq. ( |4.35| ) becomes 



Sreg(^) — S REG (0) 



PiP 



-Aq- 



p-2 
EA 



00 , cos(o;r) 
dr — ; — 



T 



(4.36) 



leading to 
t 



REG 



- S R eg(0) « p(p - l)Ag^ 2 r[l - a] sin 



«7T 



a; 



(4.37) 



where T[x] is the Gamma function. The higher order terms neglected in Eq. ( |4.35| ) 
give rise to terms vanishing as higher powers of uj. Likewise, 



5regM - ?reg(0) « 2Ar[l - a] sin 



«7T 



a; 



ia-1 



(4.38) 



\uJ\ 



Replacing Eqs. (|437D and fli~38D in Eq. and comparing terms of order 1 ' ■ 1 " 1 

we obtain 

pip- 1 ) „v-io l -y' 



which compared with Eq. ([4.24|) gives 



x Vl 



p-2 



y' 



(4.39) 



(4.40) 



the marginal value of x p . We thus see that the marginally stable state is the only 
one compatible with a gapless spectrum. Conversely, the equilibrium state must 
necessarily have a gap in its excitation spectrum. 

In order to determine the values of the exponent a and the amplitude A, we 
proceed by keeping only the first term on the rhs of Eq. (|4.35| ), 

Sreg(^) - Sr E g(0) « P ^ P ~ l \ ll 2 (ftusoW - <?reg(0)) , (4.41) 

and replace it in Eq. ( |4.23| ) noting that (?reg(0) = 5^(0) — j3q EA = j3mq EA /(p — 2) (cf. 
Eqs. ([4T8D and The result is 



^1 P(p 



p-l _ P[P 1 l ,,!>-2;, 

V p-2 A """ <BA 2 



-/3mg| 



feEGl^fcJ 



1 . 



(4.42) 



The solution of this quadratic equation is 



<?REG(^fcJ 



^ + ^Pmq^-\u Jk \T-y\/f + 



2p(p-l) 
P-2 



P-i 

EA 



p(p - l)q EA 2 



(4.43) 
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This result is only valid in the limit ujk — > where it reduces to 
taM ~ ?d(0) - /3<? EA 



(&(0) - /3g EA ) 1 / 2 g E r p/2 v /27R^Tj] 



\Uk\ 



(AAA) 



Thus, the exponent giving the asymptotic decay of g<j(r) in Eq. (|4.34! ) is a = 2. 
Substituting Eqs. ( |4.19| ), (|4.20|) and ( |4.4Q ) in Eq. ( |4.44|) and analytically continuing 
to real frequencies we obtain the exact result 



lim 

w— >0 U) 



1 

7f 



p(p - 1) 



3/4 



(4.45) 



One can easily show that the amplitude A is obtained by multiplying the rhs of the 
equation above by n. 

It is interesting to notice that a linear excitation spectrum has also been found 
in the case of the SU(AT) Heisenberg SG model [32]. However, in our model this 
gapless spectrum is not a consequence of Goldstone's theorem as the model does 
not possess any continuous symmetry. 

4.4.5 An approximate solution of the dynamical equations 

The results of the previous section suggest a simple approximate solution of the com- 
plete set of equations that becomes exact in the low-frequency and low-temperature 
limit. Once qd{oj) is known, the equation of state, q EA (T,T), must be determined 
from the normalization condition 



1 



X >)coth(/3u;/2) . (4.46) 

JO 71 



If we assume that the integral is dominated by the low frequencies, we can compute 
it by replacing x"^) by the analytic continuation of Eq. ( 4.43|) . The result at T = 
is 

A VT 

' &a = —7 ^174 > ( 4 - 47 ) 



3tt 



p(p - 2 /2 



which reduces to 



r 



3tt 
T 



pip - l)gL" 2 /2 



1/2 



Qe 



(4.48) 



This equation has the same form as the corresponding one obtained with the static 
approximation, Eq. (|3.33|); only the coefficient is different. For p = 3, we find that 



the maximum coupling for the marginal SG state, that is also the critical coupling at 
zero temperature, is T d (0) = 2.75 compared to r rf (0)=3.97 in the static approxima- 
tion and the estimate 1^(0) = 3.1 from the numerical calculations to be described 
below. It is easy to check that, at finite but low temperature, the corrections to 
Eq. flP7D are of 0{T 2 ). Therefore, 



g EA (T,T) =g EA (0,r) 1-0{T 



(4.49) 
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as T -> 0. 

We shall compare quantitatively the low-frequency approximation with the nu- 
merical solution of the exact equations in Section ||. 



4.5 Thermodynamic functions 



The free-energy per spin may be calculated by substituting our Ansatz for the order- 
parameter matrix in Eqs. ( [2.23| ) and (|2.24j) . After some algebra we obtain: 



PS 



-|log[l-(l- m)y'\ 
P rfi 



m 



rn 



^-(m - l)g| A - y / drq p d (T) - + In 



log- 

(3z 



i-y' 



(1 — m)y' 



2 sinh 



5> 



UJ 



+ z q d (u) 



+ 5? 



UJ 

T 



+ z q d {uj) 



Tzjl 
1 



(4.50) 



The entropy, internal energy and specific heat are given by 

d(Pf) 



b ~ dT 



U 



dp 



_dU 
v dT 



(4.51) 



respectively. 

The derivatives needed to compute S and U must be taken only with respect 
to the explicit T-dependence of the free-energy. Contributions from terms of the 
form d(Pf)/dq EA x dq EA /dT and alike vanish because of the stationarity condition. 
A convenient way to proceed consists in observing that q d (r) is a dimensionless 
function and, as such, it can only depend on dimensionless variables. We can choose 
them as follows 

q d (r, J,r,T) = q d (r/P,Pz,Plr/J)- 

It follows that 



ff 3 



£ Cds f d (s,Pz,pj,T/J) 



4 



(4.52) 
(4.53) 



and its explicit derivative with respect to P is 
dU 2 _ 

dp " 2 ./,, 

Similarly, 



f f 1 dsq p d (s,Pz,pj,T/J) 



1 

2 Jo 



drq p d (r) 



drq d (r)e 



tU! k T 



P I dsq d (s)e^ ks = Pq d (k) 



(4.54) 



(4.55) 



The explicit (3 dependence of q d {uj) is thus a multiplicative factor. Finally, we 
observe that duj^jdp = —u k /p. Using these considerations, we derived the following 
expression for the internal energy: 



+ z 



(4.56) 



We analyze separately the PM and SG phases. 
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4.5.1 Internal energy of the paramagnetic phase 

The equation of motion can be recast as 



1 - 



+ z q d {u) k ) = -E(u k )q d (u k ) 



(4.57) 



replacing this expression in the last term of Eq. (|4.56 ) and using 

P 



1 



^2T,(u k )q d (uj ri 



" drnr)q d {r) = V -fdr q p (r) 
o l Jo 



(4.58) 



we obtain a very simple expression for the paramagnetic internal energy: 

z p + 2 fP 



dr q p d (r) . 



(4.59) 



Since we have shown that a gapless spectrum is only possible in the marginally 
stable SG state, there must be a gap E G in the spectrum of the paramagnetic phase. 
Therefore the specific heat in the PM phase vanishes exponentially at low tempera- 
ture: 

C v oc exp(-£ G /T) . (4.60) 
4.5.2 Internal energy in the spin-glass phase 

The same procedure can be applied to the SG phase, using this time the equation 
of motion in the form 



<jJt 



where 



/ P = l + 



q d (u k ) = l + (f p - l)<L fc ,o 



X , 



m(l + Xp) m 2 (l + Xp) 
By substituting Eq. ( [4.61|) into Eq. ( [4.56D we obtain 



^sg = | - 13m q p EA - J dr (q p (r) - q EA ) 



where (3m is given by Eq. (|4. 1 9|) and 



z = -Pmq^ +Xp + S REG (^ = 0) 



x , , 



and 



t REG (uJ = 0) = | J dr (q p d \t) - q^ A 



(4.61) 



(4.62) 



(4.63) 



(4.64) 



(4.65) 



In the equilibrium SG state the specific heat vanishes exponentially as its exci- 
tation spectrum has a gap. 
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Since / is not an extremum with respect to m in the marginal state, it would seem 
that the term d(/3f)/dm x dm/dT should appear in the calculation of its internal 
energy. However, it can be shown that in the correct expression of the energy 
density of this state, that coincides with the one obtained in the real-time dynamics 
calculation, such contributions are absent. Hence, we shall use Eq. ( |4.56 ) also in the 
marginal case. Since now q REG (T) = qdij) — q^A ~ t 1 at T = 0, we expect C v to 
have a power-law dependence upon T. Indeed, since at low temperature the integral 
in Eq. ( [4.63| ) is dominated by the long-time behavior of ^(t), we can evaluate 
the latter by Fourier transforming Eq. ( 4.43|) . The resulting expression depends 
on T only through q EA . It follows from Eq. ( |4.49| ) that, in the marginal SG state, 
C v oc T. This power-law behavior is reminiscent of the temperature dependence of 
low-temperature glasses that has been intensively studied experimentally [O . 



5 Numerical solution 

In this Section we present numerical solutions obtained using an algorithm that 
solves the equations of Section |j. All the results reported in this Section where 
obtained for p = 3. 

5.1 The paramagnetic phase 



The PM phase is described by Eq. ( |4.2|) subject to the spherical constraint (|2.27|) . 



At each temperature, the equations are solved iteratively as a function of T, starting 
from high and low values of this parameter. In these limits we can find analytical 
perturbative solutions that are used as starting points for the iteration procedure. 
For each value of T, the input is the converged solution obtained for the previous 
value. 



For (3 < (3 P ~ 6, Eq. (O) has only one solution that is connected continuously 
to both the high and low-r limits. Fig. ^|(a) shows the function ^(r) for (3 = 4 
and several values of T. The correlation function in imaginary time decreases mono- 
tonically as a function of r in the interval [0, (3/2]. The initial slope increases with 
T. For large T the decay is exponential, reflecting the presence of a well developed 
gap in the excitation spectrum. Fig. ||(b) represents the T-dependence of the free- 
energy per spin for the same value of (3. For T 3> 1, it approaches a T-independent 
asymptote, / = T/8. 

For (3 > (3 P , the solutions obtained for T ^> 1 and r C 1 are still unique but 
they are not continuously connected. The solution that derives from that for r ^> 1 
(PMi) can only be followed down to a critical coupling T = T c2 , where it ceases to 
exist. Conversely, the solution that exists for Y <C 1 (pm 2 ) can be followed only up 
to r = r d- Since T c2 < T cl the two solutions coexist in the interval T c2 < T < r cl . 
Figure |7| shows the two types of solutions for (3 = 8 and several values of T. For 
this temperature, the two paramagnets coexist in the region 2.8 < T < 3.6. The 
corresponding free-energies are shown in Fig. (b). The free-energies cross at a 
point intermediate between T c2 and r c i. Notice that, below the crossing point, the 
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Figure 6: (a) The function q d {r) for the PM solution at (3 = 4 and T = 1,2,3,4,5 
from top to bottom, (b) The PM free-energy as a function of T at (3 = 4. The dashed 
line is the asymptotic value / = T/8. 



classical PM has the lowest free-energy. These results are reminiscent of those found 
within the static approximation with PMi and PM 2 corresponding to the solutions 
with q d = and qd = q(( of Section |3], respectively. 

The two PM states are qualitatively different. In the quantum PM, ^(r) decreases 
exponentially in < r < (3/2. In the classical PM, after a transient, qd{i~) levels- 
off to a r- independent value that decreases with increasing T. An analysis of the 
connection between qdij) and x"( u ) given in Eq. ( |2.30| ) allows us to understand 
these differences in terms of the spectral properties of the two states. We find that 
the exponential decay of the correlation function of the quantum paramagnet reflects 
the presence of a gap in ~x"{uj) and that the plateau that appears in the case of the 
classical paramagnet reveals the presence of a central peak with a narrow width 
Au; <C T in the gap. The hight of the plateau measures the fraction of the total 
spectral weight under the peak. 



5.2 The spin-glass phase 

The equations that describe the SG phase contain an extra parameter, the break 
point m. The solutions are found as follows. Fixing (3 and m, we compute q EA , (jd{0) 
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Figure 7: (a) qd{j) for the PM solutions at /3 = 8. Solid lines: classical PM solutions. 
T = 1,2, 3, 3.5 from top to bottom. Dashed lines: quantum PM solutions. V = 
2.8,3,3.5,4 from top to bottom, (b) Solid Lines: free-energies of PM solutions as 
functions of T at — 8. The solid circles are the points at T c i = 3.6 and T C 2 = 2.8 
where the solutions cease to exist. Dashed line: the asymptotic value / = T/8. 

and z' from Eqs. (|4.19| ), (|4 . 2 0|) and ( }4.24|) , respectively. The value of x p is chosen 
according to whether we want to study the equilibrium or the marginal SG states. 
We then solve Eq. ( |4.23|) iteratively varying T, starting from T = 0, until a value is 
found that satisfies the spherical constraint. For each temperature lower than the 
classical transition temperature, this procedure is repeated for different values of m. 

This procedure allows us to determine a function m = m(T,T). As shown in 
Fig. | for (3 = 20, m(T, T) has two branches that meet at T = r max (T). Physical 
values of m lie on the branch that satisfies that m(T, T = 0) = ^class 

(T), the 

classical break point. There are no solutions of Eq. (|4 . 2 3|) for T > r max (T) which 
may be identified as the value of the coupling above which quantum fluctuations 
destroy the SG phase. 

We found that, as in the static approximation, there exists a temperature T* 
such that, for /?>/?*, m^^iT) < 1. In all cases q EA is finite at r max (T). This fact 
and Eq. (|4.19|) imply that lim-r-^o ^max(^) = 0, meaning that replica symmetry is 
restored at the quantum critical point. The same feature was previously found in a 
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Figure 8: The breaking point m as a function of V for p = 3 and (3 = 20, both 
for the equilibrium (solid line) and marginally stable (dashed line) SG state. For 
the equilibrium state (r MAX , m MAX ) ~ (3.21,0.7). For the marginally stable SG state 

(rMAX) m MAx) — (3.12,0.5). 



study of the SU(A/") Heisenberg model [|32jl . 

Fig. H shows the regular part of the auto-correlation function (cf. Eq. Q4.21| )), 
for p = 4 and (3 = 12 and values of T in the SG phase. It may be seen that q REG ( r ) 
decays more rapidly in the equilibrium state than in the marginally stable state. 
Analysis of the curves of Fig. [5] for (3 = 12 shows that, for 1 <C r <C (3/2, q REG (r) 
decays exponentially in the first case but q REG (r) oc T ~ 2 i n the second one. This 
follows again from the differences in the excitation spectra of the two states. 

In the case of the marginally stable solution, q REG (^~) ma y also De computed 
using the low-frequency approximation of Section [4.4.5| . The exact and approximate 
results for (3 = 20 are shown in Fig. ID]. It can be seen that the agreement between 
the two is very good. 



5.3 Equilibrium phase diagram 

In this section we discuss the phase diagram that results from the analysis of the 
equilibrium solutions. The case of the marginally stable solutions will be discussed 
in the next Section. As within the static approximation, all the PM states that 
we found are locally stable. Their free-energies must be compared with that of 
the SG state in order to construct the equilibrium phase diagram. As before, there 
exists a temperature interval T* < T < T P within which there are two first order 
transitions as a function of T, one from the quantum to the classical PM, followed 
at a lower coupling by the SG transition. However, for p = 3, the case for which we 
have performed detailed numerical calculations, T P — T* is very small and, in the 
following, we only consider in detail the cases T > T P and T < T*. Figures |TT|(a) 
and (b) show the T-dependence of the PM and SG free-energies for p = 3 computed 
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Figure 9: Regular part of the auto-correlation function in the SG phase. Full and 
dashed lines represent the equilibrium and marginally stable states, respectively. 
Upper curves: (3 — 4, Y = 2.35. Lower curves: (5 = 12, V — 1.75. 
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Figure 10: Comparison of the exact numerical solution for q R EG{ T ) (solid lines) and 
the approximation of Section |4.4.5| (symbols) for (3 = 20. Upper curves: T = 2.0. 
Lower curves T = 3.0. 
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using Eq. ( 2.24 ) for two temperatures, above T P and below T*. Solid lines and 
symbols represent the PM and SG solutions, respectively. The curves end at the 
point where the corresponding solution disappears. It may be seen that for T > T P 
the free-energies of the two states merge precisely at T C (T) where m—1. Below the 
critical point, / SG > / PM meaning that the SG solution maximizes the free-energy. 
As in the static approximation, stability arguments do not exclude the PM solution 
below T c as a metastable state but we may resort to the arguments given in Section 
^ to argue that the continuation of the PM solution into the SG phase is unphysical 
also in this case. 

Below T* we must compare the free-energies of the two PM solutions with that 
of the SG state to choose the most favorable phase. In this temperature range, the 
classical paramagnet has the same unphysical properties that we found in Section ^| 
in the static approximation: the ground-state energy and susceptibility diverge and 
its free energy does not intersect that of the ordered state. Therefore, we discard 
it in favor of PMi even if its free-energy is higher. The free-energies of the SG and 
PMx states cross at T c < r max as shown in the inset in Fig. [TT](b) . As in the static 
approximation, / SG < / PM for T <T* whereas / SG > / PM for T > T*. The SG and PMi 
solutions extend beyond the point where their free-energies cross and there is phase 
coexistence. Now, both q EA and m are discontinuous at r c . The SG transition is thus 
first order with latent heat and discontinuous susceptibility (see below). The phase 
diagram resulting from this analysis is represented in Fig. [E^ (thin lines). The flat 
section is the first-order line, which does not exhibit the spurious reentrant behavior 
found within the static approximation (Fig. ||). In addition, in contrast with the 
approximate phase diagram, the transition line has zero-slope at T = 0, which is a 
consequence of the third law of thermodynamics as shown in Ref. [|l| 



We also computed the T-dependence of q EA and the static susceptibility, 



X 



/ dr [q d (r) - (1 - m)q EA ] , (5.1) 
Jo 



as functions of T for the p = 3 model. The results are displayed in Fig. |13| The 
susceptibility has a cusp at T c for T > T* and a discontinuity for T < T*. The 
dotted lines correspond to the regions of metastability. Notice that the susceptibility 
is higher on the SG side of the transition, an unusual result. This follows from the 
fact that the gap in the excitation spectrum of the PM state is wider than that of the 
SG state. For the same reason the entropy of the PM is lower leading to a negative 
latent heat at the transition. 



As shown in the lower pannel of Fig. 13, the order parameter decreases rapidly 



with increasing T and is reduced by a factor of two half-way from the transition, 
showing that quantum fluctuations are quite strong in this system. 

As mentioned above, early experiments jlj on LiHo a: Yi_ :1 .F4 in a transverse mag- 
netic field suggested that the second order SG transition observed for moderate fields 
could become first order for H m H c (0), the field above which the system remains 
PM down to zero temperature. The recent observation of hysteresis effects in the 
transverse field-dependence of the susceptibility give further support to this idea M. 
While the model that we study here is not intended to describe microscopically this 
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Figure 11: Free-energies of the different PM (solid lines) and SG (symbols) phases 
above T P (a) and below T* (b). The inset in panel (b) shows in detail the crossing 
of the free-energies at the critical point for T < T*. 




Figure 12: Static (thin lines) and dynamic (thick lines) phase diagrams of the p-spin 
model for p = 3. Solid and dashed lines represent second and first order transitions, 
respectively. 
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compound, it captures some of its phenomenology. However, important differences 
between the predictions of this model and experiment exist, the most important 
ones being the fact that, in the model, where the transition is second order, it is 
not accompanied by a divergence of X3 an d that the jump of x across the first order 
transition line has a sign opposite to that observed experimentally. 

6 Connections between the marginally stable state 
and real time dynamics 

As is well known, the dynamics of classical spin-glasses becomes non-stationary be- 
low a dynamic transition temperature Td that may or may not coincide with T g . In 
the low-temperature phase, the systems shows aging, i.e., the correlation function 
C(t,t') depends both on the time-difference (t - if) and on if [ffOj pil], |42| . Aging 
has experimental manifestations such as a dependence of the ac-susceptibility on 
frequency and the time elapsed after a quench. Theoretical studies of aging effects 
in quantum systems have recently appeared in the literature pUL 03, EH HH]. Some 



systems where quantum fluctuations might play a role, like magnetic nanoparti- 



cles |47| and disordered 2D-electronic systems j|8j , were investigated experimentally. 
A search for these effects in the Li(Ho/Y)F4 system is under way ||. 

The existence of connections between dynamic and static calculations in classi- 
cal disordered systems is known since the work of Sompolinsky on the equilibrium 



dynamics of the SK model [49|. Apart from some conceptual problems that have 
been extensively discussed in the literature pOfl , Sompolinsky noted that his dy- 
namic solution for the classical SK model was marginally stable. Kirkpatrick and 



Thirumalai BT[ showed that, in the classical mean-field p-spin-glass model, Td can 
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be calculated in a fully static replica approach if the break point m is computed 
by demanding that the spin-glass susceptibility Xsg diverges, instead of asking that 
the free-energy be an extremum with respect to variations of this parameter. In 
technical terms, this prescription corresponds to setting the replicon eigenvalue to 
zero |52| and goes under the name of condition of marginal stability. It has already 
been used within the Matsubara approach to quantum disordered problems on the 
basis of the physical assumptions of having a gapless solution |53|, |32| . 

In this Section we compare the results obtained above using the marginal stability 
criterion with the outcome of a full dynamical calculation of the same system weakly 
coupled to a quantum bath made of an ensemble of harmonic oscillators ]2Uj . The 
Hamiltonian of the coupled system is 



H = H SYST + H BATli + aH mT , (6.1) 



with -ff S YST given in (^2). H BATK is the Hamiltonian of an ensemble of quantum har- 
monic oscillators with an Ohmic spectral density. H mT represents a linear coupling 
between the spins and the oscillators' displacements and it is controlled by a cou- 
pling constant a. The oscillators are in equilibrium at temperature T at all times 
and act as a reservoir. 

The real-time evolution of the system was derived using the Schwinger-Keldysh 
formalism and it is described by coupled differential equations for C(t, t') and R(t, t'), 
the symmetrized correlation and response functions of the system, respectively |2(J . 
In the glassy phase, C exhibits a two-time dependence of the form: 

C(t + t w , t w ) = C ST (t) + C AG (t + t w , t w ) . (6.2) 

The first term gives the time-dependent evolution for time- differences t <C t w and 
varies between 1 — q EA to 0. In this time regime, the second term is constant and 
equal to q EA but decays from q EA to in a waiting-time dependent manner when 
t^> t w . More formally, the Edwards- Anderson parameter is defined as 

q EA = lim lim C(t + t w ,t w ) . (6.3) 

t — >oo t w — >oo 



The response function can also be written by separating two terms as in ( |6.2|) 

R(t + t w , t w ) = R ST (t) + R AG (t + t w , t w ) . (6.4) 

The first term is invariant under time-translations and it is related to the first term 
in (|6.2j ) by the quantum fluctuation-dissipation theorem (FDT). The second term 
is related to the second term in ( |6.2| ) by a modified quantum FDT relation where 
the temperature of the bath T is replaced by an effective temperature T EFF = T/X 
which defines X, the FDT violation factor In the long-i w limit, 



P (+ \+ + \ X dC AG (t + t w ,t') 



T dt' 



(6.5) 



t'—t 



For the model under consideration, weakly coupled to a bath, the FDT violation 
factor and q EA are given by |2(J 
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I3X = ^^=0)^-2)^,/, (6.6) 
1 = ^f^C- 2 -R?> = 0). (6.7) 



These equations have the same form as Eqs. ( |4.32|) and (|4.33| ) with m and <? RE g(0) 



replaced by X and R st (uj = 0), respectively. It will be shown below that, in the 
limit a — > 0, the equation obeyed by R ST (u) reduces to the analytic continuation 



of Eq. Q4.231 ). Therefore, in this limit, R st (lu = 0) = (/reg (^fc = 0). It follows 
from this fact and Eqs. (^32|)-(^33|) and CT-(p77|) that m = X. Although the 
coincidence between the values of X and m for the marginal SG state has been 
noticed several times for classical models, this is the first explicit evidence of its 
validity in a quantum problem. 

This correspondence allows us to identify the dynamic transition line T^T) as 
the boundary in the T — T plane where the marginally stable solution exists. T^T) 
is a piecewise curve formed by a line on which X = 1 that joins (T^, 0) and (T*, T*) 
and a line on which X < 1 that joins (T*,T*) and (0, T c ). The dynamic phase 
diagram for p = 3 is shown in Fig. [12] (thick lines). It is qualitatively similar to 
that obtained in the static approximation. In particular, m is discontinuous across 
the dashed line and Td lies always above T c . Notice that the two transition lines 
come very close to one another for T ~ T* and seem to touch at T*. For T < T*, 



m varies continuously along T&(T) and vanishes at the quantum critical point [55 
This has a consequence of potential interest for experiment: at low temperatures, 
FDT violations are predicted to appear suddenly as r d is crossed coming from the 
high T region. 

The equivalence between the replica theory supplemented with the marginality 
condition (rtmc) and the dynamical approach goes beyond the equality of m and 
X. It also holds for the expectation values of operators calculated within rtmc and 
the long-time limit of the same observables in the coupled dynamical system, for 
a —>■ 0. It must be stressed that the weak-coupling limit must be taken after the 
thermodynamic and long-time limits [[54]] . More precisely, one can establish that for 



an operator O depending only on the coordinates and momenta of the system, 

(0) RTMC = lim lim lim {0(t)) DYN , (6.8) 

a— >0 t— too N^oo 

where ( • • • ) rtmc is the outcome of the rtmc calculation. 

To illustrate this point, we consider Eq. fl4.23|) fixing p = 3 for simplicity. We 
define the real-time quantities 

^RTMcM = gHBo(iWfc) Ifc, _ MJ+i0 f , #RTMC(£) = IZo ^ ^RTMC (w) eXp ( ~lUt) , 



(6.9) 



Srtmc(^) = T, REG (iu k ) i(J ^ +i0+ ; s rtmc(£) = f-oo ^ £ RTM0 (u;) exp(-iut) 
and introduce the Lehmann representations 



- , . /°° c(uImi? RTMC (w) ~ [°° dulmE KTMC (uj) 
<?R EG (^fc) = / : , T, RBG (iu k ) = / ; . (6.10) 

-00 7T LO — lUJk J -00 7T LO — lUlk 
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It follows from Eqs. ( p.9| ) and ( |6.10| ) that 

^RTMc(t) = 0(t) JZo % ImR RTUC (Lu) sm(cjt) , 
SrecW = 0(t) JZo V ImS RTMC (^) sin(^) . 



(6.11) 



Analytically continuing Eq. (|4.23 ) and Fourier transforming the resulting equation 
we arrive at 



1 d 2 
f dp 



+ z n 



-Rrtmc(^) — + / dt E RTMC (t — t )-RrtMc(^ ) 
Jo 



(6.12) 



where the limits on the integral on the right-hand side follow from Eq. (|6.11 ) and 

^rtmc = z' + S REG(0) = R~^ mc (uj = 0) + J dt' S RTMC (t'). (6.13) 

The last equality follows from the combined use of Eqs. ( 4.24 ), Eqs. ( |4.19| ) and 
( |4.33|) . Rewriting Eq. ( [4.25 ) for p = 3 in the form 



2(3 



^2qwa(i^>m) Qwcii^k - i^m) + 3g EA g R E G (^fc) 



(6.14) 



and using Eq. QS.lUp , we find after some algebra 

, , 3 r°° du' 



Im S T 



Crtmc(^ ) im -Rrtmc(^ — 

-oo 71 



(6.15) 



where we have introduced the Fourier transform of the symmetrized correlation 
function 



Crtmc(^) = Im R rtmc (uj) coth + 2n q EA 5(u). 



It follows that 

^rtmc(^) = 3 0{t)C RTMC (t) -Rrtmc(^) = — 3 Im 



CrTMC (t) 2 "^RTMC (t) 



(6.16) 



(6.17) 



Equations ( |6.12| ), ( |6.13| ) and ( |6.17| ) are equivalent to Eqs. (6.6), (6.16), (7.13) and 
(3.27) of Ref. [ 20 1 specialized to p = 3 and taken in the limit a — > 0. It follows that 

-RrtmcW = lim lim R DYN (t + t w ,t w ) = limi? ST (t) . (6.18) 

a^0iu,-+oo a— >0 

The equality of the rtmc response function and the stationary part of the dynamical 
response function implies that all the quantities that can be expressed in terms of 
them are also equal. In particular, it can be shown that the asymptotic dynamic 
energy- density of the system in contact with the bath in the weak coupling limit is 
precisely given by Eq. (|4.63 ). 
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7 Conclusions and discussion 



In this paper we presented a detailed study of the properties of the quantum spherical 
p-spin-glass model in thermodynamic equilibrium and in the marginally stable state. 

We solved the equations that describe the different phases of the system nu- 
merically and also using various approximation schemes that give valuable physical 
insights, in particular, on the delicate issue of handling the multiplicity of solutions 
of the equations. 

We established explicitly a connection between the states obtained in the quan- 
tum replica calculation with the marginality condition and those that result form 
the non-equilibrium dynamics of the system weakly coupled to an environment in 
thermal equilibrium. 

We determined the phase diagram of the model and showed that, in all cases, 
quantum fluctuations drive the SG transition first-order at sufficiently low temper- 
atures. A tricritical point separates the two transition lines. This shows that the 
properties of SG systems in the quantum regime can be qualitatively different from 
those in the classical limit. The same feature has been first observed in the p —>■ oo 
limit of the p-spin model in a transverse field jTTJ and in a quantum model with 
multiple interactions studied in Ref. [|L2] . These studies, as well as the TAP anal- 



ysis of Ref. suggest that this phenomenon may be rather generic in quantum 
extensions of models having a discontinuous transition in the classical limit. 

We notice, however, that Ritort studied the phase transition in the random or- 
thogonal model (ROM) |56| in a transverse field within the static approximation and 



found no evidence of the existence of a tricritical point [p7 |. In this case, an expan 



sion of the free-energy in powers of m — 1 leads to a vanishing Edwards- Anderson 



parameter at the quantum critical point (0, T c ) |j58[ . It would be interesting to inves- 
tigate whether this behavior subsists beyond the static approximation, i. e. whether 
a tricritical point at finite temperature appears in the exact solution. 

A quantum particle (or a manifold) in an infinite dimensional random potential 
has a discontinuous classical transition when the potential has short-range corre- 
lations [62]. Therefore, the quantum extensions analyzed by Goldschmidt and 



Giamarchi and Le Doussal [53] are hence candidates to exhibit a first order transition 
in the limit of zero temperature. 

So far, studies of quantized models with continuous classical transitions have 
found continuous transitions close to the quantum critical point. However, we believe 
that this issue should be explored further. Full solutions of the low-temperature 
phase of these models do not yet exist. They would need a full replica symmetry 
breaking Ansatz for the non-diagonal terms in the replica matrix. In our opinion the 
soft spin quantum model whose real-time dynamics has been studied by Kennett 
and Chamon |46j is the simplest choice to investigate this question. 

Many of the results of this paper can also be obtained using the TAP ap- 
proach [ 23H - This formalism, that avoids the use of replicas, allows us to understand 
the organization of the metastable states and to interpret the dynamic behavior 
of glassy systems in terms of them. The relationship between the quantum TAP 
equations and the Matsubara approach for this model was analyzed In Ref. [BB[ . A 
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relationship between (q m , m, qdi^)) in the replica approach and (g EQ , S, C(r)) in the 
TAP approach was obtained where S is the energy density of a TAP solution and 
C(r) the imaginary-time dependent correlation function. Different choices for m, 
such as the equilibrium condition the marginality condition or others, correspond 
to different choices for £: the equilibrium value, the value reached dynamically or 
some intermediate one, respectively. 

In Section [T2] we mentioned that this model can be viewed as a description of 
a closed polymer embedded in an infinite dimensional space with a random poten- 
tial |25|, [27]]. The line tension favors a linear configuration of the polymer while 
the disorder makes the polymer wander and search for configurations that are fa- 
vorable from an energetic point of view. As seen above, besides a transition line 
separating a liquid-like phase from a glassy-like phase, we expect a transition line 
between two different liquid phases corresponding to the two paramagnetic solutions 
that we found in the spin-glass model. One of the liquid phases is characterized by 
a large value of the correlation of two monomers separated at maximum distance 
L/2, with L the total length of the polymer. (In the static approximation we called 
this solution .) In the polymer language, the other solution corresponds to a state 
in which the correlation of two monomers at distance L/2 is small and close to zero 
at high values of T. These two configurations can be interpreted as representing a 
coiled polymer (q^) and a linear-like polymer (q^). These results seem reasonable 
since the solution q^ appears at high values of T, i.e. when the polymer is very 
flexible, while the solution q^ appears at low values of V when the polymer is more 
rigid. The nature of the transition from the liquid-like to the glassy-like phase is dif- 
ferent in the two situations. This transition is expected to occur on a sizable region 
of phase space only for p ^> 1 which corresponds to short range correlations of the 
random potential. It would be interesting to check whether such a crossover occurs 
in finite dimensional models of random directed polymers and in experiments. 

This work can be extended in several directions. One of them is the question of 
how the choice of particular initial conditions affect the real-time dynamics of these 
disordered quantum models. Up to now, only the case of random initial conditions, 
which simulate an instantaneous quench from the high-temperature phase, was dis- 
cussed in the literature [F21J. An interesting problem is the relaxation dynamics 
within metastable states that are solutions of the quantum TAP equations [p3||, in 



the manner of Refs. [5(| and [p9|| . To solve this problem one has to represent the 



generating functional of the correlation functions in the form of a path integral on 
a time contour that mixes real and imaginary times. The derivation of the relevant 
equations and their solution are not straightforward and need, as input, the solution 
of the imaginary-time equilibrium equations studied in this paper. We shall discuss 
this problem elsewhere [|60|| . Another open problem is how the coupling to a strong 
quantum environment modifies the statics and dynamics of such disordered models, 
an issue that we shall discuss in a separate publication |61 . 
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Appendix 

In this appendix, we demonstrate that, for the model studied in this paper, the 
extremization equations for a 2-step RSB Ansatz collapse into the ones for the 1-step 
RSB. This procedure can be iterated for a k-step RSB in order to prove that the 
1-step RSB scheme is exact. 

Within a /c-step RSB Ansatz, the n replicas are organized into njm\ "1-families" 
of mi elements each. The m\ replicas in a "1-family" are then further organized into 
mi/m2 "2-families" of m 2 elements each, and so on until the level of "A;-families" is 
reached. We use the convention that all replicas belong to the same "0-family" so 
that m = n and rrik+i = 1. Off diagonal elements of a A;-step RSB matrix Q a b for 
which a and b belong to the same "/-family" are labeled q t . A 2-step RSB matrix 
Q is then characterized (at zero magnetic field, that is, q = 0) by two off-diagonal 
elements, q\ and q2, and two break points, m\ and m 2 , such that, in the n — > limit, 
< mi < m 2 < 1. The only terms in the free-energy density that depend on these 
parameters are given by 



1 



log(g d (0)//3 - ((g))) + — log ( 1 + mi - <h 



I m2 - 1 log ^ ^( )/^~^ 2 



m 2 *\qM/P-((q)), 



m, & V q d (0)/P-{(q))J 
- j(m 2 - l)ql - ^(mi - m 2 )q{ , (A.l) 



where ((g)) = g 2 + m 2 (gi —52)- Hence, the extremization equations in order to q±, 
g 2 , mi and m 2 are given by 

<h + 1-<, ! ; 1 = , (A.2) 



[g~d(0) + (m 2 - l)(3q 2 + ( mi - m 2 )/3g 1 ](g~ d (0) - (3{(q))) 2 



1 - 



Pqi Pqi - 9d(0) 

(9d(0) + (m 2 - l)Pq 2 + {mi - m 2 )pq 1 ) q d {0) - j3q 2 

P 



X g d (0)-/3((g)) + l Arl 



0, 



(A.3) 
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9 

mf 



log 



&(O)-0« ff » 



k ?d(0) - (1 - m 2 )[3q 2 + (mi - m 2 )f3qi 
1 



+ 



and 



ml 



log 



»™i ?d(0) - (1 - m 2 )/3q 2 + (mi - m 2 )/?gi 
Qi-Q2 ( 



(A.4) 



+ 











(A.5) 



g d (0)-/3((g» V^(0)-/5((g}} + mi/3 gi m 2 J 

respectively. As can be easily seen, for q± = q 2 and mi = m 2 (or m 2 = 1), the above 
set of equations yield the 1-step result. In order to demonstrate that the 1-step 
solution is unique, we now compute the linear combinations — px Eq. ( |A.4|) +/? 2 (/iX 
Eq. (|A~2|) and (3 2 q 1 x Eq. ( |Q|) -q 2 x Eq. Q) -px Eq. flA3[) and define x p = 
(/5?i)/?d(0) and y = (/3q2)/qd(0), to obtain the following equations: 

-,,2 



+ 



(1 + (m 2 - l)y p + (mi - m 2 )x p )(l - y p - m 2 (x p - y p )) 
1 px p 



mi 1 + (m 2 - l)y p + (m x - m 2 ) 



(A.6) 



p 



P_ log / 1 + (m 2 - l)y p + (mi - m 2 )x p N 



and 



.1 



y p -m 2 (x p -y p ) 
P 



y p - m 2 (x p - y,,) 
x p (p + 1) 



+ - 



m 



log 



1 + (m 2 - l)y p + (mi - m 2 )x p 
'\-y p - m 2 (x p - y p ) x 



+ 



y P , _p_ 
- y P "^2 



i - y P 



o. 



respectively. At this point, it is useful to define the parameters 

m 2 (x p - y p ) 



Up -- 
and 

■i 

and rewrite eqs. (|A.7|) and ( |A.7|) as 

„,2 



i -y P 

miXp 

i -y P 



(l-u p + V p )(l - Up) 



+ 



VpP 



Up + v p 



p log 1 + 



and 



'm 2 Wp(l +p) m 2 \ «p 

mi 1 — Up + v p nil I 1 — u % 



plog 



u.,. 



0. 



(A.7) 

(A.8) 
(A.9) 

(A.10) 
(A.ll) 



respectively. The only real solution to the above equations is u p = and v p equals 
its 1-step value (for example, v p = 1.81696 for p — 3). Therefore, we can conclude 
that the 2-step RSB scheme yields back the 1-step one. 
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